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I. INTRODUCTION 



Bessel functions of the first kind (J n ) are among the most ubiquitous special functions 
of mathematical physics. Their properties are described in several monographs^-, in the 
encyclopedic Abramowitz & Stegun 5 as well as in every advanced mathematics textbook. 

The generating function identity 



£ UZ) W U (1) 



n=— oo 



is the key formula used in discussing the frequency modulation (FM) of the electric field of 
a laser. In fact, a pure sinusoidal FM reads as 

E(t) = 5° ( e ^t+Msin(fH)) + c c ) ^ (2) 

where ul is the laser carrier frequency, Q is the modulation frequency, M the modulation 
index and MQ the modulation depth. Using (TjQ) with w = e lQt one finds the well-known 
sidebands structure 

E(t) = ^ Yl J "( M ) e i( " i+nn) * + c.c. (3) 



2 

ra=— 00 

In some applications, for instance studying the response of a two-level system to a FM field, 
one has to deal with expressions like 

_ ^ J n (M)J n _ s (M) 

n=—oo 

where 7 represents a resonance linewidth. Even if the sum could be done exactly using the 
Newberger sum rule&, the result is quite involved (the details are given in Section DTLA]). 

Alternatively, in the physical interesting limit e = < 1 (see below for more details), 

a geometric expansion gives 



I - —//-(- I n 2 + 0(e 3 ) 



(5) 



A s =- V J n (M)J n „ s (M) l-^n-(-) 

reducing the problem to the evaluation of sums of the form 

+00 

B Ks = n k J n (M)J n „ s (M) (6) 



for integer values of k > and s. To the authors' best knowledge these sums are not 
addressed in the literature in the general case although special values are known, for instance 
B 0fi = 1 from 9.1.76 of ] . 

Finding a closed expression for B k jS is the main goal of the present work. The paper 
is organized as follows. In Section [TT] the main results are derived and discussed, while in 
Section [TTT1 some applications are pointed out. Finally Section HVl contains the conclusions. 

II. RESULTS 

We start the derivation of the main result by rewriting ([ID in a slightly different form 

oo 

f{y,0) = * iv ** 9 = E J n (y)e* n9 . (7) 

n=—oo 

Deriving the left-hand side k times, one gets 

f ik \y,0) = ^^= Pk (y,9)f(y,9) (8) 
where Pk(y, 0) is a polynomial in cos# and sin# and can be easily recast in the form 

oo 

Pk (y,8)= E C k , n (y)e ine , (9) 

n=— oo 

with Ck, n {y) polynomials in y. An explicit form for these coefficients is worked out in 
Appendix [A] Combining everything one finds 

oo / oo \ 

f (k) (y,o)= E E c Uy)Jn-M) ^ (10) 

n=— oo \q=—oo / 

On the other hand, deriving the right-hand side of ((?]) one obtains 

oo 

f^(y,6)= E ^n) k J n (y)e ind . (11) 

n=— oo 

An explicit comparison of (TTTj) with ( TTUj) is performed in Appendix [B] Now consider the 
expression F(yi,y 2 ,8) = f^(yi,9)f(y 2 ,0). A first form can be obtained using OH]) and ([7]) 
as follows 

F( Vl ,y 2 ,9)=p k (y 1 ,9) e J ^ + ^ sin9 



E c k , n (yi)e ine \ E J m(yi + y2)e l 

=— oo / \m=— oo 

D / OO \ 

E E C k , q - m (y 1 )J m (y 1 + y 2 ) )e^ e . 



(12) 



q=— oo \m=— oo 



I III t 



(13) 

qt) 



In the meantime an equivalent formula can be written using (jTTj) 

too \ / oo 

Y {in) k J n { yi )J n6 \ J ™(^)e 
n=— oo / \m=— oo 

oo / oo \ 

= Y ( ( tn } k Jn ^ J <2-«(^) ) e 

q=—oo \n=— oo / 

and, given the completeness of the set {e tq6 , q = 0, ±1, ±2, . . .}, the coefficients in both 
expansions must be equal 

oo oo 
Y C k,q-m(yi) Jm{Vl+V2) = (i 7i) fe J n (?/i) J g _ n (y 2 ) . (14) 

m=— oo n=— oo 

This is the central result of the paper and can be considered to be a generalization of the 
Jn addition formula 



formula 9.1.75], which, in fact, is obtained for k = noticing that 
Co,n = 9 (see Appendix [A}. Besides, with a proper choice of y\ and y 2 , different sum rules 
can be obtained. In fact, substituting y 1 = —y 2 = y in f[T4"l) and remembering that J„(0) = 
one finds 

oo 

C k ,q(y)= Y ( in "> k Jn ^ J i-ni.-y) 

n=~ oo 
oo 

= Y (in) k J n (y)(-l) q - n J q -n(y) (15) 

n=— oo 

oo 

= Y ( in "> k Jn ^ Jn -v(y)- 

n=— oo 

This gives an answer for the initial problem, namely 

Bk,s = -j: Ck,s (16) 

and the Bk, s for the first k values can be read from Appendix IA1 

The companion case y± — y 2 — y is worked out with similar algebraic elaborations leading 
to the formula 

Y nk J M J n-M = —T- E C k ^ m {y)J m {2y). (17) 

n=— oo m=— oo 

A. Generalizations 

Many generalizations are possible once one dives into the world of generalized Bessel 
functions 4 . Let us discuss some cases. 
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Consider the function 

oo 

e i( x cos e +y sme) = J n (x, y) e* n \ (18) 

n=— oo 

which arises in the context of general sinusoidal modulation. The J n (x,y) are a kind of 
"generalized" Bessel functions, which can be given in terms of the usual Bessel functions by 
means of simple algebraic manipulations (see also equation (11.10) in Dattoli 4 ) 

oo 

J n {x,y)= i q J q (x) Jn- q (y). (19) 

q=— oo 

The same steps as before lead to the sum rule 

oo 

2 ^2 nJ n (x,y) J n -g(x,y) = (y + ix) 8 q>1 + (y - ix) 8 q ^ x . (20) 

n=— oo 

Furthermore, /c-order derivatives lead to sum rules with n k which the interested reader can 
easily work out. 

The second generalization is achieved considering 

oo 

el ( yiS m0 +s/2 sin2fl) _ £ J n { yi ,y 2 )e ind . (21) 

n=— oo 

These J n {y\-,V2) "generalized" Bessel functions can be expressed as 

oo 

Jn(yi,y2)= J l(y^ J n-2q(yi), (22) 

q= — oo 

and are real functions for real values of yi (see chapter 2 of Dattoli 4 - for further details). By 
means of the same device one gets 

oo 

^ n J n (yi,y 2 )J n -s(yi, 2/2) = 7^(^,1 + <^,-i) + 2/2(^,2 + S a - 2 )- (23) 

n=— 00 

Also in this case higher derivatives lead to other sum rules. 

In a complete general setting, one could consider the problem of a FM signal with no 
special properties, i.e. 



c 



i<f>(t) 



J2 G n e inQt (24) 



where the only requirement is that (f)(t) is a real and periodic well-behaved function (a more 
precise statement is given below) 

00 

<j>(t) = <t>ne inm . (25) 



The G n coefficients are related in a very complicated way to multiple sums of products of 
ordinary Bessel functions, but the explicit form is not important here. They satisfy the 
"conservation of energy" sum rule 

oo 

Gn G* n _ s = 5 S , Q . (26) 

n=— oo 

Taking the time derivative of ( 124)) and doing the same steps as above we get 

oo 

nG n G* n _ s = isct> s . (27) 

n=— oo 

Equation f[2T|) holds provided that <fi(t) can be expanded as a meaningful Fourier series, i.e. 

oo 

0(t) = inQ0 n e innt (28) 

n=— oo 

must be convergent. 

As before, higher order derivatives will give higher order sum rules. 

III. APPLICATIONS 

An example is considered where formula (fl~6j) comes in handy. Let us model an optical 
resonance as a damped harmonic oscillator of unit mass, forced by the electric field of the 
light beam 

z + yz + u%z= fe +iu) \ (29) 

where, as usual, the physical oscillator displacement is x(t) = 3l(z(t)) and all other symbols 
have a clear meaning. Such a simple and analytical model avoids all the complications tied 
to the solution and lets us concentrate on the sum rule. In fact, after the transient, the 
solution settles into the well-known 

z(t) = 2 J 2M . , (30) 

and the work done by the light field per time unit is 

dw f\ iut Ms ( iu}eiut -iue~ iuJt \ 

- V( e + e ) , . 2 73 I ~ . + ,. 2 ,,o t- • (31) 



dt 4 \loq — tu 2 + ijoj co 2 — to 2 — i r yuj / 

The experimental apparatus usually averages out high frequencies and thus the measured 
signal is 

{ dt } ~ 2^\u 2 -u 2 + i 1 u)~2 (u 2 - u 2 ) 2 + w 2 7 2 ' ( ' 



which, for nearly resonant light uj = coq + 5, 5 <C coq, becomes 



[ dt ' 2 



1 1 1 

+ 



A 3 



O 



(33) 



_ 7 1 + A 2 2(1 + A 2 ) 2 coo V^5/j 

where A = 25/7. It is easily seen that the leading term shows the well-known Lorentzian 
form. 

Replacing the forcing term with a frequency-modulated one 



i ui t y i (ui t+M sin Qt) 



J n(M)e l 



0J n t 



(34) 



on the right-hand side of fl29|) we get the exact result 

f 



dW 

v dt 1 



■3 



yv * y> u n J n (M) J n _ s (M) | ^,0,, 



-00 \n=— 00 



£jg - W 2 + i 7Wr , 



(35) 



where o; n = c^o + & + w fi. To proceed any further it is necessary to simplify the expression 
F = tv n / (ojq — uj\ + i^ujn). In the optical range it is not restrictive to assume u ^> 7, 5 thus 
getting 



26 + 2 nVt - i 7 



1 + 



(5 + nfi) 5 



25 + 2 nfi -17 w 



Considering only the leading term in F we are faced in f[35|) with the sums 

J„(Af)J„_,(M) 



1 E - 

ry ^ 1 



27 <— ' i + iA + 2nifi/ 7 

n=— 00 ' 



(36) 



(37) 



which are very similar to (j3J). To proceed any further with the geometric expansion of the 
denominator, we notice that the product J n (M) J n _ s (M) shows the well-known "bridge" 
structure and it is exponentially small when \n\ > Nmax ~ 2 M. So if 

2N MAX tt/ 1 



|2AWfi/ 7 | < 1 
and the sum can be elaborated as follows 

Nmax 



1 -HA 



< 1 



(38) 



1 E- 

i-l ^ 1 



J n (M)J n . s (M) 



1 -r - 

Z7 ^— ' 1 



J n {M)J n _ s {M) 



27 ' l + iA + 2nif2/7 27 z — ' 1 + 2A + 2 

n=-oo n=-N MAX 



* N M a: 

1 E 

27 _ 



J n {M)J n _ s {M) 



n=-N MAX 



1 + iA 



[l + ne + nV + 0(e 3 )] (39) 



00 

1 E 

«7 ^— ' 



J n (M)J n . s (M) 
1 + iA 



[l + ne + nV + 0(e' 



where e = — i (2fl/j) /(1 + iA) and the last row is justified because of the exponentially small 
nature of the terms included. 

The inner sum in (135]) can now worked out by means of the sum rules developed above, 
and, after some straightforward algebra, one finds 



where the terms in parentheses, besides the d.c. component, represent the higher harmonic 
modulated absorption. Notice how the A dependence of the (2MQ/ / y) n term is related to 
the real and imaginary part of n-th derivative of 1/(1 + zA). It seems worth stressing that 
the validity of (jUJ) is limited to modulation depths smaller than the linewidth as pointed 
out in (1581) . Usually in typical experimental conditions one is interested in the first harmonic 
signal, catching it, for instance, by means of a lock-in amplifier. The above formula shows 
that it is more convenient to use the biggest value of the modulation depth M Q compatible 
with (1381). 

More generally, using an arbitrary but limited modulation signal as in (1251) one finds 



ulation depth" a = \ujmax — u min \ is put in evidence to compare with previous formulae. 

A. Comparison with the Newberger sum rule 

As stated in the Introduction, the quantity A s can be summed by means of the Newberger- 
sum rule. The result is 




(40) 




(41) 



where in the last row the instantaneous frequency u(t) = <fi(t) is introduced and the "mod- 



A s 



(-l) s 7r 7 /fi 



J s - i7 /n(M) J 47/ n(M) 



s > 



(42) 



7 sinh(7T7/f2) 



S 



where one has to deal with complex-order Bessel functions. Even if this is an exact result, 
the physics is obscured by the complexity of the formula. For instance, it is not easy to 
extract the behaviour when rj = fi/7 is small. Instead using (jSJ) with s = 1 we quickly found 

Ai = ~ i V ~ ^ V 2 + ^ (1 + 3 M 2 /4) i V 3 + 0(r/ 4 ) (43) 
27 27 27 

To proceed further using (142]) the product of the Bessel functions must be developed 5 
(formula 9.1.14), and, making use of the properties of the Euler Gamma function^ (formula 
6.1.31), one gets 

(-ir /MVA (-^+a)! /- 1 * 1 

which, to the authors' best knowledge, is an original elaboration. This formula is still 
complicated, but better suited to extracting asymptotic behaviour. In fact for s — 1 we 
have 

M f 1 3M 2 1 1 \ 

1 ~ ~2^ [l-i-y/tl ~ 4 2-^/fi T + TtT^F + " 7 (45) 

= ~ iV ~ V 2 + + 3 M 2 /4) i ^ 3 + O(^). 
27 27 27 

IV. CONCLUSIONS 

We have derived some new classes of sum rules and recursion relations obeyed by Bessel 
functions of the first kind. The results have broad applications in the physics of modulated 
excitation for instance in the case of light-matter interaction. A comparison of our results 
with an elaborated form of the Newberger sum rule is also given. 



Appendix A: Determination of Ck >n 

The starting point is the recursion obeyed by the pk polynomials. In fact, it is easily seen 
that po = 1, and assuming (JHJ) then follows 

Pk+i(y, 0) = (JIq + i y cos 9J p k (y, 9). (Al) 
Projecting this equation as shown in (j9j) one obtains 

y 

Cfc+i,n = i n Ck jTl + i — (Ck jn+ i + Ck,n-i) with Co, n = 8n,o- (A2) 



9 



To the authors' best knowledge this recursion cannot be solved in terms of known special 
functions, even if the homogeneous right-hand side, namely i n Ck, n + i § (Cfc,n+i + Ck,n-i) = 
0, is satisfied by the J n (y) itself. The first few loops give 

^ 1 , n = |(+l + -l) 

4^2,n = V - (20 + +2 + -2) + V - (+1 - -.1) 
^3,n = £ (+3 + -3) + ^ (+2 - -2) + 



4^n = ^(+4 + -4) + ^(+3--3)4 
z 4 lo 4 



^ + !|!)(+2 + -2) + 
J/ 

4 2 
3y 4 ?/ 2 \ 

The recursion is straightforward, but tedious. A simplification is achieved thinking of ( IA2j) 
as a matrix vector product and implementing it in any computer algebra system. 

Alternatively, an approach based on the Faa di Bruno formula can be developed. The 
Faa di Bruno formula generalizes the chain rule to higher derivatives and can be stated as 
follows 

d k 



dx 



t(miA \-m k )t „/\\ TT / I* V / 



e -r-r — f / (mi+ - +mfc) ( y (x))n f- 

f— ' mi! m 2 ! • • • m k \ W 

{mj} 3=1 x 



where the sum is over all /c-tuples of non- negative integers mj, . . . , m& satisfying the con- 
straint 1 mi + 2 m 2 + 3 m 3 + • • • + & = k. 
Applying that device to our problem we find 

d k 

= -iy sine i y sin 6 

Pk d6 k 

u\ k ( A4 ) 

= Y tt , ., {iy) m Y[(sm^e) m ^ 
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where m = X]j>i m i- Next using shv 2 - 7 l ' 9 = (— 1)^ 1 - ) cos 6 and sin*- 2 ^ 6 = (— 1)^ sin# we 
get 

_ u 

P* = E TT ^(^ri-l^sin^cos^ (A5) 

where = Ej>o( m 2+4j + ™2+4i+i), a = m 2j and 6 = m 2i _i = m - a. Finally 

expressing the trigonometric functions as 

sin a £ CO s b 6 = -^J2J2(- i y( a ) ( b )e^ m - 2r ^ e (A6) 

r=0 r'=0 \ / \ / 

and substituting in (1A5|) we find 

{m 3 } X1 J J J r=0 \ / \ 2 / 

with m > n > —m, where we follow the convention that the binomial coefficients are zero if 
the lower index is not an integer or negative or larger than the upper one. 



Appendix B: Some new recursion relations 

By comparing the two different forms (TlO|) and (TIT]) of the same quantity, a relation 
between the coefficients results 

q k J q (y) = nC k M J q - n (y). (Bl) 

n=—oo 

The k = case is the identity J q (y) = J q {y)- For k = 1 we obtain the well- known- (formula 
9.1.27 first row) three term recursion relation 

2qJ q (y) = y [J q+1 (y) + J q -i(y)]- (B2) 

The interesting and, to the authors' best knowledge, new relations are found for k > 2. 
In fact for k = 2 one gets 

£ - ■/,(*) 

7 2 (B3) 

y y 

= g [ J 9+i(z/) - Jff-i(2/)] - "J [^+2(2/) + ^-2(2/)]- 
For higher values of k, other relations are easily obtained using the coefficients of Appendix 

m 
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